Structural basis of reactivation of oncogenic p53 mutants by a small molecule: methylene quinuclidinone (MQ).

In response to genotoxic stress, the tumor suppressor p53 acts as a transcription factor by regulating the expression of genes critical for cancer prevention. Mutations in the gene encoding p53 are associated with cancer development. PRIMA-1 and eprenetapopt (APR-246/PRIMA-1MET) are small molecules that are converted into the biologically active compound, methylene quinuclidinone (MQ), shown to reactivate mutant p53 by binding covalently to cysteine residues. Here, we investigate the structural basis of mutant p53 reactivation by MQ based on a series of high-resolution crystal structures of cancer-related and wild-type p53 core domains bound to MQ in their free state and in complexes with their DNA response elements. Our data demonstrate that MQ binds to several cysteine residues located at the surface of the core domain. The structures reveal a large diversity in MQ interaction modes that stabilize p53 and its complexes with DNA, leading to a common global effect that is pertinent to the restoration of non-functional p53 proteins.

Expression of a recombinant, 4'-Phosphopantetheinylated, active M. tuberculosis fatty acid synthase I in E. coli.

BACKGROUND: Fatty acid synthase 1 (FAS I) from Mycobacterium tuberculosis (Mtb) is an essential protein and a promising drug target. FAS I is a multi-functional, multi-domain protein that is organized as a large (1.9 MDa) homohexameric complex. Acyl intermediates produced during fatty acid elongation are attached covalently to an acyl carrier protein (ACP) domain. This domain is activated by the transfer of a 4'-Phosphopantetheine (4'-PP, also termed P-pant) group from CoA to ACP catalyzed by a 4'-PP transferase, termed acyl carrier protein synthase (AcpS). METHODS: In order to obtain an activated FAS I in E. coli, we transformed E. coli with tagged Mtb fas1 and acpS genes encoded by a separate plasmid. We induced the expression of Mtb FAS I following induction of AcpS expression. FAS I was purified by Strep-Tactin affinity chromatography. RESULTS: Activation of Mtb FAS I was confirmed by the identification of a bound P-pant group on serine at position 1808 by mass spectrometry. The purified FAS I displayed biochemical activity shown by spectrophotometric analysis of NADPH oxidation and by CoA production, using the Ellman reaction. The purified Mtb FAS I forms a hexameric complex shown by negative staining and cryo-EM. CONCLUSION: Purified hexameric and active Mtb FAS I is required for binding and drug inhibition studies and for structure-function analysis of this enzyme. This relatively simple and short procedure for Mtb FAS I production should facilitate studies of this enzyme.

Structure of Type-I Mycobacterium tuberculosis fatty acid synthase at 3.3 Å resolution.

Tuberculosis (TB) is a devastating and rapidly spreading disease caused by Mycobacterium tuberculosis (Mtb). Therapy requires prolonged treatment with a combination of multiple agents and interruptions in the treatment regimen result in emergence and spread of multi-drug resistant (MDR) Mtb strains. MDR Mtb poses a significant global health problem, calling for urgent development of novel drugs to combat TB. Here, we report the 3.3 Å resolution structure of the ~2 MDa type-I fatty acid synthase (FAS-I) from Mtb, determined by single particle cryo-EM. Mtb FAS-I is an essential enzymatic complex that contributes to the virulence of Mtb, and thus a prime target for anti-TB drugs. The structural information for Mtb FAS-I we have obtained enables computer-based drug discovery approaches, and the resolution achieved by cryo-EM is sufficient for elucidating inhibition mechanisms by putative small molecular weight inhibitors.

New Insights into the Role of DNA Shape on Its Recognition by p53 Proteins.

The tumor suppressor p53 acts as a transcription factor recognizing diverse DNA response elements (REs). Previous structural studies of p53-DNA complexes revealed non-canonical Hoogsteen geometry of A/T base pairs at conserved CATG motifs leading to changes in DNA shape and its interface with p53. To study the effects of DNA shape on binding characteristics, we designed REs with modified base pairs "locked" into either Hoogsteen or Watson-Crick form. Here we present crystal structures of these complexes and their thermodynamic and kinetic parameters, demonstrating that complexes with Hoogsteen base pairs are stabilized relative to those with all-Watson-Crick base pairs. CATG motifs are abundant in p53REs such as GADD45 and p53R2 related to cell-cycle arrest and DNA repair. The high-resolution structures of these complexes validate their propensity to adopt the unique Hoogsteen-induced structure, thus providing insights into the functional role of DNA shape and broadening the mechanisms that contribute to DNA recognition by proteins.

Post-translational regulation of p53 function through 20S proteasome-mediated cleavage.

The tumor suppressor p53 is a transcription factor that regulates the expression of a range of target genes in response to cellular stress. Adding to the complexity of understanding its cellular function is that in addition to the full-length protein, several p53 isoforms are produced in humans, harboring diverse expression patterns and functionalities. One isoform, Δ40p53, which lacks the first transactivation domain including the binding region for the negative regulator MDM2, was shown to be a product of alternative translation initiation. Here we report the discovery of an alternative cellular mechanism for Δ40p53 formation. We show that the 20S proteasome specifically cleaves the full-length protein (FLp53) to generate the Δ40p53 isoform. Moreover, we demonstrate that a dimer of FLp53 interacts with a Δ40p53 dimer, creating a functional hetero-tetramer. Consequently, the co-expression of both isoforms attenuates the transcriptional activity of FLp53 in a dominant negative manner. Finally, we demonstrate that following oxidative stress, at the time when the 20S proteasome becomes the major degradation machinery and FLp53 is activated, the formation of Δ40p53 is enhanced, creating a negative feedback loop that balances FLp53 activation. Overall, our results suggest that Δ40p53 can be generated by a 20S proteasome-mediated post-translational mechanism so as to control p53 function. More generally, the discovery of a specific cleavage function for the 20S proteasome may represent a more general cellular regulatory mechanism to produce proteins with distinct functional properties.

Diverse p53/DNA binding modes expand the repertoire of p53 response elements.

The tumor suppressor protein p53 acts as a transcription factor, binding sequence-specifically to defined DNA sites, thereby activating the expression of genes leading to diverse cellular outcomes. Canonical p53 response elements (REs) are made of two decameric half-sites separated by a variable number of base pairs (spacers). Fifty percent of all validated p53 REs contain spacers between 1 and 18 bp; however, their functional significance is unclear at present. Here, we show that p53 forms two different tetrameric complexes with consensus or natural REs, both with long spacers: a fully specific complex where two p53 dimers bind to two specific half-sites, and a hemispecific complex where one dimer binds to a specific half-site and the second binds to an adjacent spacer sequence. The two types of complexes have comparable binding affinity and specificity, as judged from binding competition against bulk genomic DNA. Structural analysis of the p53 REs in solution shows that these sites are not bent in both their free and p53-bound states when the two half-sites are either abutting or separated by spacers. Cell-based assay supports the physiological relevance of our findings. We propose that p53 REs with long spacers comprise separate specific half-sites that can lead to several different tetrameric complexes. This finding expands the universe of p53 binding sites and demonstrates that even isolated p53 half-sites can form tetrameric complexes. Moreover, it explains the manner in which p53 binds to clusters of more than one canonical binding site, common in many natural REs.

AcpM, the meromycolate extension acyl carrier protein of Mycobacterium tuberculosis, is activated by the 4'-phosphopantetheinyl transferase PptT, a potential target of the multistep mycolic acid biosynthesis.

Modification of acyl carrier proteins (ACP) or domains by the covalent binding of a 4'-phosphopantetheine (4'-PP) moiety is a fundamental condition for activation of fatty acid synthases (FASes) and polyketide synthases (PKSes). Binding of 4'-PP is mediated by 4' phosphopantetheinyl transfersases (PPTases). Mycobacterium tuberculosis (Mtb) possesses two essential PPTases: acyl carrier protein synthase (Mtb AcpS), which activates the multidomain fatty acid synthase I (FAS I), and Mtb PptT, an Sfp-type broad spectrum PPTase that activates PKSes. To date, it has not been determined which of the two Mtb PPTases, AcpS or PptT, activates the meromycolate extension ACP, Mtb AcpM, en route to the production of mycolic acids, the main components of the mycobacterial cell wall. In this study, we tested the enzymatic activation of a highly purified Mtb apo-AcpM to Mtb holo-AcpM by either Mtb PptT or Mtb AcpS. By using SDS-PAGE band shift assay and mass spectrometry analysis, we found that Mtb PptT is the PPTase that activates Mtb AcpM. We measured the catalytic activity of Mtb PptT toward CoA, using an activation assay of a blue pigment synthase, BpsA (a nonribosomal peptide synthase, NRPS). BpsA activation by Mtb PptT was inhibited by Mtb apo-AcpM through competition for CoA, in accord with Mtb AcpM activation. A structural model of the putative interaction between Mtb PptT and Mtb AcpM suggests that both hydrophobic and electrostatic interactions stabilize this complex. To conclude, activation of Mtb AcpM by Mtb PptT reveals a potential target of the multistep mycolic acid biosynthesis.

Structural studies of p53 inactivation by DNA-contact mutations and its rescue by suppressor mutations via alternative protein-DNA interactions.

A p53 hot-spot mutation found frequently in human cancer is the replacement of R273 by histidine or cysteine residues resulting in p53 loss of function as a tumor suppressor. These mutants can be reactivated by the incorporation of second-site suppressor mutations. Here, we present high-resolution crystal structures of the p53 core domains of the cancer-related proteins, the rescued proteins and their complexes with DNA. The structures show that inactivation of p53 results from the incapacity of the mutated residues to form stabilizing interactions with the DNA backbone, and that reactivation is achieved through alternative interactions formed by the suppressor mutations. Detailed structural and computational analysis demonstrates that the rescued p53 complexes are not fully restored in terms of DNA structure and its interface with p53. Contrary to our previously studied wild-type (wt) p53-DNA complexes showing non-canonical Hoogsteen A/T base pairs of the DNA helix that lead to local minor-groove narrowing and enhanced electrostatic interactions with p53, the current structures display Watson-Crick base pairs associated with direct or water-mediated hydrogen bonds with p53 at the minor groove. These findings highlight the pivotal role played by R273 residues in supporting the unique geometry of the DNA target and its sequence-specific complex with p53.

Diversity in DNA recognition by p53 revealed by crystal structures with Hoogsteen base pairs.

p53 binds as a tetramer to DNA targets consisting of two decameric half-sites separated by a variable spacer. Here we present high-resolution crystal structures of complexes between p53 core-domain tetramers and DNA targets consisting of contiguous half-sites. In contrast to previously reported p53-DNA complexes that show standard Watson-Crick base pairs, the newly reported structures show noncanonical Hoogsteen base-pairing geometry at the central A-T doublet of each half-site. Structural and computational analyses show that the Hoogsteen geometry distinctly modulates the B-DNA helix in terms of local shape and electrostatic potential, which, together with the contiguous DNA configuration, results in enhanced protein-DNA and protein-protein interactions compared to noncontiguous half-sites. Our results suggest a mechanism relating spacer length to protein-DNA binding affinity. Our findings also expand the current understanding of protein-DNA recognition and establish the structural and chemical properties of Hoogsteen base pairs as the basis for a novel mode of sequence readout.

Anomalous DNA binding by E2 regulatory protein driven by spacer sequence TATA.

We have investigated the anomalously weak binding of human papillomavirus (HPV) regulatory protein E2 to a DNA target containing the spacer sequence TATA. Experiments in magnesium (Mg(2+)) and calcium (Ca(2+)) ion buffers revealed a marked reduction in cutting by DNase I at the CpG sequence in the protein-binding site 3' to the TATA spacer sequence, Studies of the cation dependence of DNA-E2 affinities showed that upon E2 binding the TATA sequence releases approximately twice as many Mg(2+) ions as the average of the other spacer sequences. Binding experiments for TATA spacer relative to ATAT showed that in potassium ion (K(+)) the E2 affinity of the two sequences is nearly equal, but the relative dissociation constant (K(d)) for TATA increases in the order K(+ )< Na(+ )< Ca(2+ )< Mg(2+). Except for Mg(2+), K(d) for TATA relative to ATAT is independent of ion concentration, whereas for Mg(2+) the affinity for TATA drops sharply as ion concentration increases. Thus, ions of increasing positive charge density increasingly distort the E2 binding site, weakening the affinity for protein. In the case of Mg(2+), additional ions are bound to TATA that require displacement for protein binding. We suggest that the TATA sequence may bias the DNA structure towards a conformation that binds the protein relatively weakly.

Structure-function analysis of the acyl carrier protein synthase (AcpS) from Mycobacterium tuberculosis.

We have solved the crystal structure of the acyl carrier protein synthase (AcpS) from Mycobacterium tuberculosis (Mtb) at 1.95 A resolution. AcpS, a 4-phosphopantetheinyl transferase, activates two distinct acyl carrier proteins (ACPs) that are present in fatty acid synthase (FAS) systems FAS-I and FAS-II, the ACP-I domain and the mycobacterial ACP-II protein (ACPM), respectively. Mtb, the causal agent of tuberculosis (TB), and all other members of the Corynebacterineae family are unique in possessing both FAS systems to produce and to elongate fatty acids to mycolic acids, the hallmark of mycobacterial cell wall. Various steps in this process are prime targets for first-line anti-TB agents. A comparison of the Mtb AcpS structure determined here with those of other AcpS proteins revealed unique structural features in Mtb AcpS, namely, the presence of an elongated helix followed by a flexible loop and a moderately electronegative surface unlike the positive surface common to other AcpSs. A structure-based sequence comparison between AcpS and its ACP substrates from various species demonstrated that the proteins of the Corynebacterineae family display high sequence conservation, forming a segregated subgroup of AcpS and ACPs. Analysis of the putative interactions between AcpS and ACPM from Mtb, based on a comparison with the complex structure from Bacillus subtilis, showed that the Mtb AcpS and ACPM lack the electrostatic complementarity observed in B. subtilis. Taken together, the common characteristic of the Corynebacterineae family is likely reflected in the participation of different residues and interactions used for binding the Mtb AcpS to ACP-I and ACPM. The distinct features and essentiality of AcpS, as well as the mode of interaction with ACPM and ACP-I in Mtb, could be exploited for the design of AcpS inhibitors, which, similarly to other inhibitors of fatty acid synthesis, are expected to be effective anti-TB-specific drugs.

Structural basis of restoring sequence-specific DNA binding and transactivation to mutant p53 by suppressor mutations.

The tumor suppressor protein p53 is mutated in more than 50% of invasive cancers. About 30% of the mutations are found in six major "hot spot" codons located in its DNA binding core domain. To gain structural insight into the deleterious effects of such mutations and their rescue by suppressor mutations, we determined the crystal structures of the p53 core domain incorporating the hot spot mutation R249S, the core domain incorporating R249S and a second-site suppressor mutation H168R (referred to as the double mutant R249S/H168R) and its sequence-specific complex with DNA and of the triple mutant R249S/H168R/T123A. The structural studies were accompanied by transactivation and apoptosis experiments. The crystal structures show that the region at the vicinity of the mutation site in the R249S mutant displays a range of conformations [wild-type (wt) and several mutant-type conformations] due to the loss of stabilizing interactions mediated by R249 in the wt protein. As a consequence, the protein surface that is critical to the formation of functional p53-DNA complexes, through protein-protein and protein-DNA interactions, is largely distorted in the mutant conformations, thus explaining the protein's "loss of function" as a transcription factor. The structure of this region is restored in both R249S/H168R and R249S/H168R/T123A and is further stabilized in the complex of R249S/H168R with DNA. Our functional data show that the introduction of H168R as a second-site suppressor mutation partially restores the transactivation capacity of the protein and that this effect is further amplified by the addition of a third-site mutation T123A. These findings together with previously reported data on wt and mutant p53 provide a structural framework for understanding p53 dysfunction as a result of oncogenic mutations and its rescue by suppressor mutations and for a potential drug design aimed at restoring wt activity to aberrant p53 proteins.

RUN-CBFbeta interaction in C. elegans: computational prediction and experimental verification.

The Runt domain proteins are eukaryotic transcription factors that regulate major developmental pathways. All members of this family contain a highly-conserved sequence-specific DNA binding domain: the Runt domain (RD). Structural and biochemical studies have shown that the Runt domain undergoes a conformational transition upon binding to DNA and that this process is regulated by an unrelated partner protein CBFbeta that enhances the DNA binding affinity of RD. Most of the reported studies on the Runt domain transcription factors were performed on proteins from mammals and Drosophila whereas very little has been known about the C. elegans RD protein, RUN, which provides the simplest model system for understanding the function of this class of transcription factors. We performed computational studies on RD domains from various species including C. elegans, Drosophila, and human, using the atom-atom contact surface area scoring method. The scoring analysis indicates that the DNA binding regulation of the C. elegans RD protein (CeRD) occurs via its interaction with a CBFbeta-like partner, as found for the human proteins, whereas a different mode of regulation may occur in the Drosophila system. Sequence, secondary structure and fold analyses of a putative CBFbeta protein identified in the C. elegans genome, CeCBFbeta, sharing a 22% identity with the human protein, predict a similar structure of this protein to that of the human CBFbeta protein. We produced the C. elegans proteins CeRD and CeCBFbeta in bacteria and confirmed their physical interaction as well as cross interactions with the corresponding human proteins. We also confirmed the structural similarity of CBFbeta and CeCBFbeta by circular dichroism analysis. The combined results suggest that a similar mechanism of regulation operates for the human and the C. elegans RD proteins despite the low sequence identity between their CBFbeta proteins and the evolutionary distance between the two systems.

Structural basis of DNA recognition by p53 tetramers.

The tumor-suppressor protein p53 is among the most effective of the cell's natural defenses against cancer. In response to cellular stress, p53 binds as a tetramer to diverse DNA targets containing two decameric half-sites, thereby activating the expression of genes involved in cell-cycle arrest or apoptosis. Here we present high-resolution crystal structures of sequence-specific complexes between the core domain of human p53 and different DNA half-sites. In all structures, four p53 molecules self-assemble on two DNA half-sites to form a tetramer that is a dimer of dimers, stabilized by protein-protein and base-stacking interactions. The protein-DNA interface varies as a function of the specific base sequence in correlation with the measured binding affinities of the complexes. The new data establish a structural framework for understanding the mechanisms of specificity, affinity, and cooperativity of DNA binding by p53 and suggest a model for its regulation by regions outside the sequence-specific DNA binding domain.

Molecular flexibility in ab initio drug docking to DNA: binding-site and binding-mode transitions in all-atom Monte Carlo simulations.

The dynamics of biological processes depend on the structure and flexibility of the interacting molecules. In particular, the conformational diversity of DNA allows for large deformations upon binding. Drug-DNA interactions are of high pharmaceutical interest since the mode of action of anticancer, antiviral, antibacterial and other drugs is directly associated with their binding to DNA. A reliable prediction of drug-DNA binding at the atomic level by molecular docking methods provides the basis for the design of new drug compounds. Here, we propose a novel Monte Carlo (MC) algorithm for drug-DNA docking that accounts for the molecular flexibility of both constituents and samples the docking geometry without any prior binding-site selection. The binding of the antimalarial drug methylene blue at the DNA minor groove with a preference of binding to AT-rich over GC-rich base sequences is obtained in MC simulations in accordance with experimental data. In addition, the transition between two drug-DNA-binding modes, intercalation and minor-groove binding, has been achieved in dependence on the DNA base sequence. The reliable ab initio prediction of drug-DNA binding achieved by our new MC docking algorithm is an important step towards a realistic description of the structure and dynamics of molecular recognition in biological systems.

Structural and energetic origins of sequence-specific DNA bending: Monte Carlo simulations of papillomavirus E2-DNA binding sites.

DNA bending is an important structural feature for indirect readout in protein-DNA recognition. The binding of papillomavirus E2 transcription factors to their DNA binding sites is associated with DNA bending, providing an attractive model system to study the origins of sequence-specific DNA bending. The consensus E2 target is of the general form ACCGN(4)CGGT with a variable four base pair region. We applied a new all-atom Monte Carlo (MC) algorithm that combines effective sampling with fast conformational equilibration. The resulting MC ensembles resemble the corresponding high-resolution crystal structures very well. Distinct bending is observed for the E2-DNA binding site with a central AATT linker in contrast to an essentially straight DNA with a central ACGT linker. Contributions of specific base pair steps to the overall bending are shown in terms of local structural parameters. The analysis of conformational substates provides new insights into the energetic origins of intrinsic DNA bending.

Structures of the DNA-binding site of Runt-domain transcription regulators.

Runt-domain (RD) proteins are transcription factors that play fundamental roles in various developmental pathways. They bind specifically to DNA sequences of the general form PyGPyGGTPy (Py = pyrimidine), through which they regulate transcription of target genes. The DNA duplex TCTGCGGTC/TGACCGCAG, incorporating the binding site for the RD transcription factors (bold), was crystallized in space group P4(3). X-ray analysis of two crystals diffracting to 1.7 and 2.0 angstroms resolution, which had slight variations in their unit-cell parameters, revealed two distinct conformations of the A-DNA helix. The two crystal structures possessed several structure and hydration features that had previously been observed in A-DNA duplexes. A comparative analysis of the present A-DNA structures and those of previously reported B-DNA crystal structures of RD-binding sites in free and protein-bound states showed the various duplexes to display several common features. Within this series, the present A-DNA duplexes adopt two conformations along the pathway from the canonical A-DNA to the B-DNA forms and the protein-bound helices display conformational features that are intermediate between those of the current A-DNA structures and that of the B-DNA-type helix of the free RD target. Based on these data and energy considerations, it is likely that the propensity of the RD-binding site to adopt the A-DNA or B-DNA conformation in solution depends on the sequence context and environmental conditions, and that the transition from either DNA form to the protein-bound conformation involves a small energy barrier.

Predicting indirect readout effects in protein-DNA interactions.

Recognition of DNA by proteins relies on direct interactions with specific DNA-functional groups, along with indirect effects that reflect variable energetics in the response of DNA sequences to twisting and bending distortions induced by proteins. Predicting indirect readout requires knowledge of the variations in DNA curvature and flexibility in the affected region, which we have determined for a series of DNA-binding sites for the E2 regulatory protein by using the cyclization kinetics method. We examined 16 sites containing different noncontacted spacer sequences, which vary by more than three orders of magnitude in binding affinity. For 15 of these sites, the variation in affinity was predicted within a factor of 3, by using experimental curvature and flexibility values and a statistical mechanical theory. The sole exception was traced to differential magnesium ion binding.

DNA recognition by the RUNX1 transcription factor is mediated by an allosteric transition in the RUNT domain and by DNA bending.

The Runt domain proteins are transcription regulators of major developmental pathways. Here we present the crystal structures of the Runt domain (RD) of the human protein RUNX1 and its DNA binding site in their free states and compare them with the published crystal structures of RD bound to DNA and to the partner protein CBFbeta. We demonstrate that (1) RD undergoes an allosteric transition upon DNA binding, which is further stabilized by CBFbeta, and that (2) the free DNA target adopts a bent-helical conformation compatible with that of the complex. These findings elucidate the mechanism by which CBFbeta enhances RD binding to DNA as well as the role of the intrinsic conformation of the DNA target in the recognition process.